Andreev Reflection in Strong Magnetic Fields 



H. Hoppe, U. Ziilicke, and Gerd Schon 

Institut fur Theoretische Festkorperphysik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 
(submitted to Phys. Rev. Lett. 26 July 1999; revised 8 November 1999) 

We have studied the interplay of Andreev reflection and cyclotron motion of quasiparticles at a 
superconductor-normal-metal interface with a strong magnetic field applied parallel to the interface. 
Bound states are formed due to the confinement introduced both by the external magnetic field and 
the superconducting gap. These bound states are a coherent superposition of electron and hole edge 
excitations similar to those realized in finite quantum-Hall samples. We find the energy spectrum 
for these Andreev edge states and calculate transport properties. 
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Rapid progress in fabrication techniques has made it 
possible to investigate phase-coherent transport in a vari- 
ety of mesoscopic conducting devices jlj. In recent years, 
the study of hybrid systems consisting of superconduc- 
tors in contact with normal metals has continued to at- 
tract considerable interest, mainly because of the novel 
effects observed in superconductor-semiconductor micro- 
junctions Many of the unusual experimental find- 
ings arise due to the phenomenon of Andreev reflection, 
i.e., the (partial) retroreflection of an electron incident on 
a superconductor (S) - normal-metal (N) interface as a 
hole |BQ . Phase coherence between the electron and hole 
states is maintained during the reflection process. Hence, 
coupled-electron-hole (Andreev) bound states || having 
energies within the superconducting gap are formed in 
mesoscopic devices with multiple interfaces, e.g., S-N-S 
systems f|, or S-N-I-N-S structures §|. (The symbol T 
denotes an insulating barrier.) Recently, measurements 
of transport across the interface between a superconduc- 
tor and a two-dimensional electron gas (2DEG) were per- 
formed with a strong magnetic field H applied in the 
direction perpendicular to the plane of the 2DEG 0. 
While the magnetic field did not exceed the critical field 
of the superconductor, it was still large enough such that 
the Landau-level quantization of the electronic motion in 
the 2DEG was important ||. With these experiments, 
a link has been established between mesoscopic super- 
conductivity and quantum-Hall physics ^ which needs 
theoretical exploration. 

In this Letter, we study a novel kind of Andreev bound 
state that is formed at a single S-N interface in a strong 
magnetic field [M. This bound state is a coherent su- 
perposition of an electron and a hole propagating along 
the interface in a new type of current-carrying edge state 
that is induced by the superconducting pair potential. 
Andreev reflection gives rise to the contribution 

2 n * 
n=l 

to the small-bias conductance, which we obtained by gen- 
eralizing the familiar Biittiker description Jllj of trans- 



port in quantum-Hall samples. In Eq. (|l|), the summation 
is over Andreev-bound-state levels that intersect with the 
chemical potential, and B n is the hole probability for a 
particular bound-state level. It turns out that n* is twice 
the number of orbital Landau levels occupied in the bulk 
of the 2DEG, and B n < 1/2 depends weakly on mag- 
netic field H for an ideal interface but oscillates strongly 
with H for a non-ideal interface. Gar can be measured 
directly as the two-terminal conductance in a S-2DEG - 
S system JtJ. Our treatment in terms of Andreev edge 
states provides a clear physical description of transport 
in such devices and explains oscillatory features in the 
conductance that were observed experimentally H and 
also obtained in previous numerical studies fll2| |. 

Let us start by recalling the classical and quantum- 
mechanical descriptions of electron dynamics in an ex- 
ternal magnetic field. When considered to be classical 
charged particles, bulk- metal electrons execute periodic 
cyclotron motion with a frequency uj c = eH/ (mc) . A 
surface that is parallel to the magnetic field interrupts 
the cyclotron orbits of nearby electrons and forces them 
to move in skipping orbits along the surface jl3) . Within 
the more adequate quantum-mechanical treatment, the 
kinetic energy for electronic motion in the plane perpen- 
dicular to the magnetic field is quantized in Landau lev- 
els |l4) which are at constant eigenvalues hui c (n+ 1/2) for 
electron states localized in the bulk but are bent upward 
in energy for states localized close to the surface jll| . Ap- 
plying the classical picture of cyclotron and skipping or- 
bits to a S-N interface (l6) , one finds that Andreev reflec- 
tion leads to electrons and holes alternating in skipping 
orbits along the interface. [See Fig. 0(a).] In what fol- 
lows, we provide a full quantum-mechanical description 
of these alternating skipping orbits in terms of current- 
carrying Andreev bound states. [See Fig. |l](b).] 

We now provide details of our calculation. A planar 
interface is considered, located at x — between a semi- 
infinite region (x < 0) occupied by a type-I superconduc- 
tor and a semi-infinite normal region (x > 0). A uniform 
magnetic field is applied in z direction, which is screened 
from the superconducting region due to the Meissner ef- 
fect. Neglecting inhomogeneities in the magnetic field 
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due to the existence of a finite penetration layer jTiJ , we 
assume an abrupt change of the magnetic-field strength 
at the interface: H(x) — HQ(x), where Q(x) is Heav- 
iside's step function. The energy spectrum of Andreev 
bound states at the S-N interface is found by solving the 
Bogoliubov-de Gennes (BdG) equation fTl], 



Ho,+ + U cxt A 

A — Ho — U ext 



= E 



, (2) 



with spatially non- uniform single-electron/hole Hamilto- 
nians H$± and pair potential JnJ A(x) = AoO(— x). 
We introduced the potential U cxt (x) — Uq S(x) to model 
scattering at the interface, and allow the effective mass 
and the Fermi energy to be different in the superconduct- 
ing and normal regions. Choosing the vector potential 
A(x) = x H Q(x) y, we have 



2mN 
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F " : • (3) 

x < 
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The operator X Py = p y £ 2 sgn(eH)/h is the guiding-center 
coordinate in x direction for cyclotron motion of elec- 
trons in the normal region, and i = *JKc/\eH\ denotes 
the magnetic length. Uniformity in the y and z directions 
suggests the separation ansatz 



u(x,y,z) = f x (x) e*» x ' 1 e lkz / \fL~T z , (4a) 
v(x, y, z) = g x (x) e^ x l e e lk * / ^L~T Z . (4b) 

The lengths L y , L z are the sample sizes in y and z di- 
rections, respectively. Solutions of Eq. (||) for the S-N 
junction are found by matching appropriate wave func- 
tions that are solutions in the normal and superconduct- 
ing regions, respectively [Q. The motion in z direction 
can trivially be accounted for by renormalizcd Fermi en- 
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h k /(2ton.s)- Non-trivial match- 
ing conditions arise only in x direction. For fixed X and 
\E\ < Aq, we have to match at x = the wave function 



fx 

fJX 



x>0 
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corresponding to a coherent superposition of an electron 
and a hole in the normal region, to that of evanescent 
excitations in the superconductor, 



with eigenvalues e± = ep 1 ^ ± E — h 2 k 2 /(2m^) and are 
assumed to be normalized to unity in the normal region. 
Hence, they are proportional to the fundamental solu- 
tions of Eq. (JsJ) that are well-behaved for x — > oo; these 

are the parabolic cylinder functions ]20| U(—j-^-, ^j^)- 
The matching conditions yield a homogeneous system of 
four linear equations for the coefficients a, b, d± whose 
secular equation determines the allowed values of E. 

It is straightforward to calculate the probability and 
charge currents for any particular Andreev-bound-state 
solution of the BdG equation (S) that is labeled by 
guiding-center coordinate X and energy E. It turns out 
that currents flow parallel to the interface. The total 
(integrated) quasiparticlc probability current is given by 



/(P) = 1 f_ dE 
x h L v dX 



(7) 



The total charge current can be written as the sum of 
three contributions, I 
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Note that I^' n ' is the current that would flow in an or- 
dinary quantum-Hall edge state ]l5[ ], i.e., due to normal 
reflection at the interface. The existence of Andreev re- 



flection is manifested in the contribution — / 



(Q,a) 



to the 



Hall current; it is proportional to the hole probability 
B(X) = J x \g x (x) | 2 . The part 7^ Q ' s) of the quasiparti- 
cle charge current is converted into a supercurrent. 

Numerical implementation of the matching procedure 
is straightforward. More detailed insight is gained, how- 
ever, when considering the limit \X\ <C \J £±/{2m-^uj 2 ) 
for which analytical progress can be made. Using an 
asymptotic form for the parabolic cylinder functions p0[ , 
the secular equation can be written as 
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The parameters 7 and A± arc defined in the usual way . 
The functions Xe± (0 solve the familiar one-dimensional 
harmonic-oscillator Schrodinger equation, 



e d 2 X e ± 

2 



de 
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Here we used the Andreev approximation j|] (E, Ao <C 
Ep ,Ej ), and the abbreviations 
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The variable v = 2 e^/(fiuj c ) coincides with the filling 
factor of quantum-Hall physics l| when the N region is 
a 2DEG. The parameter s = [6p'mN/(6p N 'ms)] 1,/2 mea- 
sures the Fermi-velocitiy mismatch for the junction, and 
w = [2mN[^o/(fi. 2 £p N ' ) )] 1 ^ 2 quantifies interface scattering. 
We discuss briefly results for two limiting cases pl[ |. 

(a) Ideal interface. In the absence of scattering at the 
S-N interface (w = 0) and for perfectly matching Fermi 
velocities (s = 1), T(a) vanishes identically. The energy 
spectrum is found from solutions of the transcendental 
equation cot (<£>+) = Ej \J Aq — E 2 . It consists of several 
bands, and states within each band are labeled by their 
guiding-center quantum number X. It turns out that 
a 2 = b 2 at any energy, and the band dispersion is 



dE 
dX 



~(N) 
•6 F 



(11) 



n i + ityj Aq — e 2 i (hu>c) 

(b) Non-ideal (S-I-N) interface at low energies. For 



E <C min{A , Hy / (2m- N X 2 )}, the dependence of tp_ 
on E and X can be neglected. We find 



E = Ar 



(2n + l)?r ± arccos(r ) - 2Xy* 2m N i i p ) /h 
q + ttA /(Huj c ) 



(12) 

where T = T{tti//2) and q = 2s/(s 2 + w 2 + 1). For 
s = 1 and w = 0, the spectrum for an ideal interface at 
small energies emerges. In the opposite limit of a very 
bad interface (s 2 + w 2 — > oo), we recover the spectrum 
of the Landau-level edge states close to a hard wall [jl5| . 
In general, Tq oscillates as a function of filling factor v. 
Hence, unlike in the ideal case, the bound-state energies 
of Eq. (|l^) vary oscillatory with v. The same is true for 
the hole probability B w b 2 < 1/2, for which we find 



B 



1 



<?7(i 



2 i + v /i- 9 V(i-r 2 ) 



(13) 



The minima in the oscillatory dependence of B on filling 
factor v occur whenever tan(7w/2) = 2u;/(l — s — w 2 ). 

Results obtained in the approximate analytical treat- 
ment sketched above are expected to be valid only as 
long as X is not too large. It turns out, however, that 
they actually provide a good description at E w for An- 
dreev levels intersecting with the Fermi energy, which are 
important for small-bias transport. In particular, we ob- 
tained a non- vanishing dispersion dE/dX close to the in- 



terface which leads to a finite Hall curr cnt4 Q ' n) -/ A . . 
[See Eqs. (||).] It is clear that, far away from the inter- 
face, i.e., for \X\ 3> r c , no coupling of electrons and holes 
via the pair potential is possible and dispersionless Lan- 
dau levels are solutions of the BdG equation. However, 
as the guiding-center coordinate X gets close to the in- 
terface, these Landau levels are bent upward and become 
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Andreev-bound-state levels for \E\ < Aq. This is seen in 
the exact numerical calculation of the spectrum (Fig. ||) , 
which also provides a crucial piece of information that is 
elusive within the approximate analytical treatment: the 
number n* of Andreev levels intersecting with the Fermi 
energy. We find that n* is twice the integer part of v/2. 

We now apply our findings to study transport in S— 
2DEG-S structures (t)]. In experiment, two S-N inter- 
faces are linked by ordinary quantum-Hall edge channels 
whose local chemical potentials differ by S/i. Generalizing 
the Buttiker formalism jll] for edge-channel transport, 
the following picture emerges. (See inset of Fig. ||.) From 
the lower edge channel, a current 51 = Sfi • e/h impinges 
on the left S-N interface. This current divides up into 
a part 6In flowing parallel to the interface in Andreev 
edge states studied above, and 5I± which flows across 
the interface. Chirality of edge states (both Andreev 
and ordinary) and conservation of quasiparticle proba- 
bility current completely determines the current parallel 
to the interface to be 5I\\ = (1 — 2B n )5[i ■ e/h. The 
upper edge channel collects 5I\\ and returns to the right 
interface, where a similar discussion applies. Hence, the 
two-terminal conductance e 5I±/5[i in the S-2DEG-S de- 
vice equals Gar [given in Eq. (El)]. Using hole ampli- 
tudes obtained from the exact numerical matching pro- 
cedure, wc determined the filling-factor dependence of 
Gar (shown in Fig. || for 2 < v < 18 and various values 
of w). As anticipated from the approximate analytical re- 
sult B n « 1/2, the ideal interface exhibits almost perfect 
conductance steps in units of 2e 2 /h. For finite scattering 
at the interface, oscillations appear in the conductance 
whose amplitude increases as the interface quality wors- 
ens. However, for certain single values of v, the ideal 
conductance is reached even at a bad interface. The lo- 
cation of minima and maxima in the field dependence 
of the conductance can be obtained from our analytical 
calculation and compare well with results of previous nu- 
merical studies Jl2][ based on a representation in terms of 
scattering states. 
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FIG. 1. Andreev bound state of an electron (solid lines) 
and a hole (dashed lines) at an S-N interface in a magnetic 
field, (a) Classical picture of electron and hole alternating in 
skipping orbits, (b) The quantum-mechanical analysis finds 
that both the electron and the hole occupy Landau-level states 
that are extended parallel to the interface. 




X (1) 

FIG. 2. Andreev-bound- state levels for the ideal S-N in- 
terface in a magnetic field, calculated numerically for v — 40 
and Ao/(Huj c ) = 2 by exactly matching solutions of the BdG 
equation for the normal and superconducting regions. 




FIG. 3. Conductance Gar, calculated according to Eq. (|lj) 
with hole probabilities B n obtained from the exact numerical 
matching procedure. We set Ao = 0.02 ?p N ' and 8 — 1. The 
value of w for each curve is indicated. The inset shows how 
the current 81 carried by the lower quantum-Hall edge chan- 
nel (solid line) is distributed at the left interface. The part 57y 
flows in Andreev edge states; it has both normal (solid line) 
and Andreev-reflection (dashed line) contributions. It is col- 
lected by the upper quantum-Hall edge channel (dot-dashed 
line) and returned to the right interface. As quasiparticle 
probability current is conserved, and backscattering prohib- 
ited due to chirality, Six equals the Andreev-reflection con- 
tribution to <57|| (dashed line) . 
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